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We consider an action that can generate fluids with three unequal stresses for 
metrics with a spacelike Killing vector. The parameters in the action are directly 
related to the stress anisotropies. The field equations following from the action 
are applied to an anisotropic cosmological expansion and an extension of the Gott- 
Hiscock cosmic string. 

I. INTRODUCTION 

General relativistic models with anisotropic stress have become increasing useful as the 



There has been considerable 



applications become more complex and precise [l,|2 

interest in anisotropic spheres [3, Is], I^, l3, 11] because of applications to stellar models 
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have generated an increased interest 



21| . Anisotropy is usually discussed 



in the context of fluid stress-energy relations, and examining the Lagrangian actions that 
generate fluid equations of state provides a simple way to view their geometric and physical 
origins. This paper discusses a simple action describing a fluid with three unequal stresses. 
The fluid is supported by a metric with a spacelike Killing symmetry. 

The vacuum Einstein fleld equations, for spacetimes with Killing vector can be gen- 
erated from an action written in terms of the Killing vector norm A = ^"^a, and twist 

= £^jafec^"V*^^. The original vacuum action was developed by Ehlers ^], Harrison |23| . 
and Geroch [24], who wrote the 3+1 Einstein equations on the 2+1 space of Killing vector 
orbits. For spacetimes with metric Qab, ( — h ++), the induced metric on the 3-manifold of 
Killing orbits is hah = Qab — ^aib/^- The action generating the vacuum fleld equations is, 
with hah replaced by 7ab := Xhab, 

Da\D''\ + DaCoD^cu^ 



2A2 



(1) 



where TZ and Da are the Ricci scalar and covariant derivative on 7 respectively, and uJa 



2 



DaOJ. This action was generalized to perfect fluids [25|, l26| by including a scalar function /C, 
and a vector function, with ^"s^ = 0. 



Sflu^d = / (i'x^in - ^-^ /C - s'^Sa) (2) 



,3 D,AD"A + DaUoD'^uo 

Using this action, a stress-energy for two isotropic equations of state can be generated. 

same action may also be applied to fluids with two 

unequal stresses. 

In this paper, which considers metrics with a spacelike Killing vector, the action Sjiuid 
is extended to completely anisotropic fluids. Previous work using this action treated /C 
and s°'Sa separately. Considering both terms simultaneously allows extensions to fluids with 
three unequal stresses. In the next section we develop the anisotropic fluid content implicit 
in S fluid, and discuss some of the equations of state and their effect on the rate-of-expansion 
of the fluid 4-velocity. Some particular metric examples are given in the third section. 



II. FLUID CONTENT 



A. Field Equations and Geometry 



The 2+1 fleld equations that follow from the action Sfiuid are 



/<-ab = \- SaSb + JClab (3j 

X^D%X-^DaX) + DaUjO^uj = (4) 

Da[^]=0 (5) 

where Da is the covariant derivative on 7. TZat = Ricci(7a6) is related to Rab = Ricci((7afe) 
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hg hi, Red = Tiab = SaSb + rCXhab (6) 

i^eRab = (7) 

KX'Rcb = (8) 



In the next section, we discuss the stress-energy related to the matter parameters, Sa and 
/C. 
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B. Stress-Energy 

A simple anisotropic fluid description contains a density and three stresses {e, Pi, i-*2, P3) 
which depend on and /C. In this section we flnd the dependence of each of the stress-energy 
components on the action parameters and show that the spatial part of is a measure of 
the stress anisotropy Pi — P2, while /C enters into the anisotropy between the (1, 3) and (2, 3) 
planes. Consider a 3+1 metric, gat, described with unit tetrad vectors [U"', e^^^, e^2)) ^"a)]) 
with spacelike Killing vector aligned with e^^y. = y/Xe^^y The general stress-energy 
tensor is assumed to be 

Tab = eUaUb + PieWei^) + P.e^^^e^^^ + P.efef^ (9) 

with the stress components carrying a tetrad index. The Ricci tensor for this stress-energy 
is 



e + Pl±P2±Pz.. ^ £ + Pi-P2-P3 (1) (1) 

Rab = 87r[ UaUb H >el 



£- Pi + P2-P3 (2)^(2) £-Pl-P2 + P3 (3)^(3) 
^ 2 ^b ^ 2 ^b r 

Using the held equation ^''^"Pab = 0, P3 is determined by the other stress-energy components 

P3 = -£ + Pi + P2 (10) 

and Ricci tensor becomes 

Rab = 87r[(Pi + P2)UaUb + {s - P2)e''^^e^b^ + {e- PM^^e^^\ (11) 

Expanding s", using Rab — SaSb + XfChab and taking scalar products, the complete stress 
energy description is 

Sa = soUa + sieW + 826^^ (12a) 

Wtts = (so)2 + (si)' + (S2)' + A/C (12b) 

IGttPi = {sof + (si)' - (S2)' - A/C (12c) 

167rP2 = (sof - (si)' + (S2)' - A/C (12d) 

167rP3 = (so)' - (si)' - (S2)' - (12e) 

S1S2 = SoSi = S0S2 = (12f) 
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where Sq, Si, and S2 are tetrad indexed. These relations between the stress-energy and 
the action parameters show how Sa and JC enter the fluid anisotropy. Using Eq. fll2cp and 
Eq. (112dll . the anisotropy in the 2+1 stress is described by the spatial components of Sa- 



87r(Pi - P2) = (sif - {32^. (13) 

The anisotropies involving the stress along the Killing direction require both Sa and /C. 

87r(Pi - P3) = sl + A/C, (14) 
87r(P2 - P3) = sl + XlC 

A particularly useful relation emerging from this formalism is the equality of the stress 
along the Killing direction with the negative Ricci scalar: IGvrPa = —P. 



C. Equations of State 



Some equations of state for the 3+1 fluid can be written down by considering values for 
Sa and /C. The conditions in Eq.( 112fl) mean that any one of Sq, Si, or S2 can be non-zero. 
In the formalism so far, the choice si or S2 equal to zero only fixes the (1, 2) index. The 
stress-energy relations for each of the choices are described in Table I (with S2 = 0) and si 
identified as the index function for stress anisotropy in the (1,2) plane. The Si = choice 
simply replaces 2 by 1 in Table 1. 



Table 1 



non-zero functions 


stress-energy relations 




s = -Pi = -P2 = -P3/3 


Sl 


e = Pi = -P2 = -P3 


So 


£ = Pi = P2 = P3 


/C, Sl 


£=-P2, P3 = 2P2 + Pi 


}C, So 


Pi = P2, P3 = -£ + 2Pi 



The fluid parameter conditions can be generally related to the value of s"Sa. For S2 = 0, 
by using Eq.( |T2l) we find 

87c{e + 2P2 - Pi) = -s'^Sa. (15) 
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If Sa is timelike, Si = 0, e + 2P2 > Pi, and Pi = P2. The timelike condition is e + Pi > 0, 
Pi = P2- If Sa is spacelike with S2 = 0, we have e + 2P2 < Pi. Since only one component of Sa 
is non zero, Sa cannot be null and we do not have e + 2P2 = Pi . Completely anisotropic fluids 
will correspond to Si (or S2) 7^ 0, /C 7^ 0. The fluid conditions for a completely anisotropic 
fluid are 

S2 = 0:e = ~P2, P3 = 2P2 + Pi (16a) 
si = 0:e = -Pi, P3 = 2Pi + P2 (16b) 



D. Energy Conditions 

There are a number of energy conditions based on the structure of the stress-energy and 
Ricci tensor In matter, the weak and dominant conditions require positive density and 
-TabU^ = eUa. With either Si or S2 non zero and a comoving observer, the strong and null 
conditions are 

Strong : RabWU^ > (17) 
-A/C=87r(Pi + P2) > (18) 
Null : TabN'^N^ > (19) 
I. N"" = U" + el^y. e + Pi> Q ^ sl> Q 
2_ = f/" + e^2) : ^ + ^2 > ^ > 
3. AT'^ = f/" + e^3) : £ + P3 > ^ (Pi + P2) > 

The weak and dominant conditions can be enforced, but the strong and some of the null 
conditions can be violated. For example, the (sq, S2 = 0) condition, e = — P2, implies 
a violation if IP2I > Pi. The violation of the strong energy condition is related to the 
development of the 4-velocity rate-of-expansion, O = VaU"", in the Raychaudhuri equation 

^ = -^0' - + U'^'cuab - RabU'^U". (20) 

dr 3 
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The action considered in this paper determines the form of the Ricci tensor. For the com- 
pletely anisotropic fluid, the Ricci tensor is 

S2 = 0:Rab = 87l[-{P^ + P^){ga, - ef )) + (Pi - P2)ei^)e[')] (21) 

= A/C(,.,-ei^)ef) + s?e«e« 
si = : P., = 87r[-(Pi + P2){gab - efef) + (P2 - P^i^^ef^] (22) 

The Raychaudhuri equation becomes 

^ = _le2 _ ^a^^fe + U'^'ujab + A/C, (23) 

dr 6 

A/C = -87r(Pi + P2). 

With zero vorticity and A/C = 0, 9 < 0, these conditions describe a decreasing rate-of- 
expansion with an initially converging or diverging time-like congruence becoming more 
focused. However, even for zero vorticity, if A/C > 0, the rate-of-change of the expansion 
can be positive, 9 > 0, and is unfocusing. 



E. Fluid Shear and Anisotropy 



The form of the stress-energy tensor, Eq.(l9]), identifies all of the spatial components as 



stress. However, equivalences 



29 



30( 1 can be used to relate the completely anisotropic form 



to fluids with shear, with the differences in the shear tensor identified as the physical cause 
of the anisotropy. The general equivalence relation for si 7^ is 



eU.U, + Pie«e« + P.e^^^e^^^ + Paef ef ^ 

e = -P2 
Pi = P- 2r7a(n) 
P2 = P - 2r7a-(22) 

Pz= P - 2^(33) 



{e + P)UaU,^Pgab-2Tia. 



ab 



(24) 
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The trace-free condition for aab provides P = (Pi + P2 + P^)/'i- Using Eq. (fT2!) the action 
parameters and K, describe the differences in the shear tensor components. 

Pi = P-2r]a(^u). (25) 
si = 167rr7[cr(22) - ^(ii)], 
A/C = 167rr7[(T(33) - a^22)]- 

Imposing an extra condition identifies an equation of state and the tetrad indexed shear 
tensor components. The evolution of spherically symmetric fluids of this type has been 
considered by Herrera et al {2], who gave a metric for a shearing expansion- free evolving 
fluid with two unequal stresses. 



III. EXAMPLES 

The stress-energy of a completely anisotropic fluid is described by 

S2 = 0:e = -P2, P3 = 2P2 + ^^1 (26a) 
s, = 0:e = -Pi, P3 = 2Pi + P2 (26b) 

P2 (Pi) must be negative for positive density. This suggests two possible application areas: 
Bianchi cosmological models, and cosmic strings with an axial tension. 



A. Anisotropic Cosmology 



Bianchi metrics have been used to model anomalies in the CMB radiation 3|, la, ll7|, ll8|, ll9| , 
with anisotropics developing due to different expansion rates along coordinate axes. A general 
metric to consider is 



-df + }?{t)dxl + c^{t)dxl + f{t)dx\ 



(27) 
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with three Hubble rates, Hi = b/b, H2 = c/c, and H3 = f /f. For a power law expansion, 
one has b = bot'^, c = CqP , f = fot^. The stress-energy is 

87re{t'^) = + (3S + -fS (28) 
87rPi(t2) = -7(7 - 1) - 6{6 - 1) - 57 
gTcP^if) = -(3{f3 - 1) - 6{6 -1)-(3S 
87rP3(t') = -1)- 7(7-1) -7/5 

The stress-energy conditions e = —P2 and P3 = 2P2 + Pi impose constraints 

7(/5 + 5)= /?(/?-!) + 5(5-1) (29) 
7/? = /?(/3-l) + 35(5-1) + (5(2/3 + 7). (30) 

Combining the constraints yields 

5(5- 1 + /? + 7) = 0. (31) 

Two cases emerge: 5 = and 5 = 1 — /3 — 7. The 5 7^ condition, when substituted back 
into Eq.( l29i) . relates /? and 7 

7 + /3 = /3' + /?7 + 7' (32) 

and results in zero energy-density. Only the 5 = case, with no expansion along the Killing 
direction, provides a fluid with non-zero density. With 5 = 0, Eq. fl29|) is 

/3(/3-l-7)=0. (33) 

When (3 = the density is zero. In order to have an anisotropic fluid with positive density, 
expansions along at least two of the axes are necessary, with b{t) = bot^ and c(t) = Cot^~^. 
For (3^0, the stress energy is 

Sne = (3{(3 - (34) 
8nPi = -iP-l)iP-2)/f 
SttPs = -(3{p - l)/t^ 
87fP, = -{(3-l){3(3-2)/f 
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This example has positive density for (3 < and /3 > 1, with (3 > 1 describing expansion. 
The rate-of-expansion for this case is 

i 0^(2^ e=(l^. (35) 

bet ^ ^ 

The congruence will unfocus for /? < 1/2. For the case of an expanding space with positive 
density the congruence is focusing. The family of radial stresses parameterized by jS is 
specially interesting. For 1 < /? < 2, the stress is a pressure, for /? = 2 it is dust, and for 
/3 > 2 it is a tension. P3 also shows this range of behavior with the P3 = dust crossover at 
(3 = 2/3, out of the physical density region. Since P3 is proportional to the negative Ricci 
scalar, the behavior of P3 also reflects a change from positive scalar curvature, through zero, 
to a negative curvature manifold. In this model, only manifolds with non zero Ricci scalar 
will have positive density. The timelike congruence will unfocus for (3 < 0. This is a positive 
density region of the parameter space with all three stresses as tensions. 

In summary (3 < unfocuses, (3 > 1 focuses. 



B. Cosmic String 



1. Metric and Stress-Energy 



Cosmic strings are of current interest both experimentally 3l|, |32|, |33|, |3J] and theoretically 



35, 



36, 



37l |. The usual cosmic string models have only a single axial tension. The simplest 



family of static 3+1 strings is described by the metric with (0, 1, 2, 3) = {t, r, z) 



-dt^ + dr^ + E^{r)dip^ + dz\ 



(36) 



The density is Sne = —E"/E, and the stresses are 



P, = P^ = 0, P, = E"/E. 



(37) 



The Gott-Hiscock(GH) static string 38|, l39| is the constant density example with a = VSvre, 
E{r) = Ci sin(ar). The matter metric matches directly to vacuum Levi-Civita with an angu- 
lar deficit and is Minkowski near the axis. There are two spacelike Killing vectors, C,(z) and 
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As an anisotropic extension to the GH string interior, consider the metric 

ds^ = -Al cos^[a(i?o - r)]dt^ + dr"^ + a"^ sin^(ar + 7)(i0^ + dz^ . (38) 

The 7 contribution to the sine argument is included for finite stress at r = 0. The stress- 
energy for this metric is 

87r£ = = -SttP,^ (39) 
SttP^ — tan[a(i?o — t)] cot(ar + 7) 
SttP^ = -a^{2 - tan[a(i?o - r)] cot(ar + 7)} 

with sin(7) ^ 0. The axial behavior of this stress-energy is especially interesting. In the 
cosmological example, several of the stress components ranged through pressure, dust, and 
tension as the metric parameter varied. In this example, the same behavior is observed 
but is linked to the radial position or, along the axis, to the size of the string. To interpret 
this example as an anisotropic cosmic string with physical tension along the axis requires 
1 < tan[ai?o] cot (7) < 2. The actual axial structure is related to the vacuum match and 
associated angular deficit. These are discussed in the next sections. 



2. Vacuum Matching 

The metric can be matched to a vacuum metric across an Israel layer aX r = Rq. The 
exterior (-I-) metric is a vacuum Levi-Civita metric with angular deficit 5 

^^^LviCivita = -dt" + dr^ + THH(t>'' + dz\ (40) 

The metric of the Israel layer is 

c^^iayer = -^Idt^ + o"^ s\\i^{aRQ + 'y)d(j)'^ + dz^. (41) 

The Levi-Civita metric is matched to the layer. The matching conditions are 

^0 = 1, a"^ sin(ai?o + 7) = ±^o5. (42) 

The stress-energy of the Israel boundary layer is calculated from jumps in the extrinsic 

curvature, Kab^ going from Levi-Civita (+) across the Israel layer to the string interior (— ), 
with the jumps calculated from < Kab > = K'^ — K~^^, K = K^. The layer stress-energy is 

- SnS^ ^<K^>-<K> g':r- (43) 



11 



The extrinsic curvatures needed to calculate the jumps are 

K+^=-6'Ro, K+ = 0, ir+=0, (44) 
= sm{aRo + 7) cos(ai?o + 7), K~ = 0, = 0. 

The surface of the anisotropic string has an Israel stress-energy content 

StcSu = -<K>, 871 S^^ = 0, SttS,, = <K> (45) 
< K > = -l/Ro - acot{aRo + 7) 

The equation of state of the GH string solution is found in the boundary layer of the 
anisotropic string. 

3. Angular Deficit 

In the GH static string, the angular deficit is related to the mass/length, /i, calculated 
from a. t = const, z = const integral of the density. For the anisotropic string, that mass is 
composed of two parts. The contribution from the string interior is 

2n Ro Ro 

f f sinfar + 7) , ,, ^ f sin (ar + 7) , 

J J a J 8tt a 

00 

= ^[cos(7) - cos(ai?o + 7)]) 
and the additional contribution from the boundary layer is 

Svra 



sin(a/?o + 7) , / d , n 
h cos(aito + 7) 



gRq 

Combining the mass densities, /i = /xi + /i2, yields 



, , sin(ai?o + 7) 
4/i = cos(7) + —-^ ^ (48) 



Substituting from the matching condition, Eq. fl42|) . the relation between the linear mass 
density and the angular deficit is 

4/i - cos (7) = ±5. (49) 
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This is not the conventional "thin string" result, 5 = 1— Ajj,. A possible explanation is 
that there is a missing potential energy associated with the shell assembly, as is found in 

. However, Futamase and Garfinkle 



a matter shell bounding Schwarzschild and vacuum 
40( 1 have pointed out that the relation between angular deficit and mass density depends 
on the matter in the string, and one could take Eg. (1491) as that relation for this anisotropic 
string. A third possibility is that there is additional, unconsidered structure at the axis. 
Noting that the metric along r = describes a 2 + 1 hypersurface, the r = axis could be 
an Israel layer boundary. This provides additional interior structure whose mass needs to be 
considered in the complete calculation. With this possibility, the actual string would be the 
interior axial structure with the anisotropic stress-energy described in this example serving 
as an atmosphere around the string. 



IV. CONCLUSIONS 



The 3+1 field equations can be written as a set of 2+1 equations on the space orthogonal 
to the Killing trajectories. In this paper we have presented a formalism which generalizes a 
simple 2+1 action for a spacelike Killing vector to describe a set of fiuids with anisotropic 
stress. The Lagrangian extension uses both a function /C, and a one-form Sadx"". The spatial 
part of Sa describes the stress anisotropy in the plane orthogonal to the Killing vector. With 
the Killing vector in the (3) direction, JC relates the energy and stress and can describe 
the anisotropies between the (1,2) and (3) planes. Sadx"" is not varied in the action. One 
advantage of using a fixed form rather than a field is its use as a modeling tool, with the 
anisotropies related to complex motions rather than new physical fields. The anisotropy is 
described in terms of stress but can be due to a number of physical mechanisms, such as 
fiuid shear. 

Two applications were considered. An anisotropic generalization of the GH cosmic string 
interior described an interior solution with positive density and anisotropic stress bounded 
by an Israel layer with the GH equation of state. The Bianchi I example described a fam- 
ily of power-law solutions containing both focusing and unfocusing expansions. The fam- 
ily of stress-energies was parameterized by a single constant, with the range of the con- 
stant describing the entire stress range: tension through dust to pressure. With the current 
strong interest in explaining anomalies in the CMB, one could consider models with a single 
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anisotropic fluid whose stress shifts from pressure, through dust to tension and focusing to 
unfocusing during a series of expansion eras. Because the stress associated with the KiUing 
symmetry is proportional to the Ricci scalar, the curvature will also evolve. In this model 
there is no expansion along the direction associated with the Killing coordinate. A com- 
pletely anisotropic expansion could be generated by considering a higher dimensional action 
and associating the Killing symmetry with a higher dimensional manifold. Both applications 
considered here underline the potential value of formally considering anisotropic stresses in 
relativity; while they form a more complicated stress-energy description, anisotropic fluids 
can have dynamic features leading to simpler matter models than those with equal stress. 
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